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According to recent progresses in the finite size scaling theory of disordered systems, thermo- 
dynamic observables are not self- averaging at critical points when the disorder is relevant in the 
Harris criterion sense. This lack of self-averageness at criticality is directly related to the distri- 
bution of pseudo-critical temperatures T c (i,L) over the ensemble of samples (i) of size L. In this 
paper, we apply this analysis to disordered Poland-Scheraga models with different loop exponents 
c, corresponding to marginal and relevant disorder. In all cases, we numerically obtain a Gaussian 
histogram of pseudo-critical temperatures T c (i, L) with mean T^ V (L) and width AT C (L). For the 
marginal case c = 1.5 corresponding to two-dimensional wetting, both the width AT C (L) and the 
shift [T c (oo) — Tc V (L)] decay as I/ -1 / 2 , so the exponent is unchanged (^random = 2 = u pure ) but 
disorder is relevant and leads to non self-averaging at criticality. For relevant disorder c = 1.75, 
the width AT C (L) and the shift [T c (oo) — T£ V {L)\ decay with the same new exponent L~ 1/Vrandom 
(where v ran dom ~ 2.7 > 2 > v vure ) and there is again no self-averaging at criticality. Finally for 
the value c = 2.15, of interest in the context of DNA denaturation, the transition is first-order in 
the pure case. In the presence of disorder, the width AT C (L) ~ L~ 1//2 dominates over the shift 
[T c (oo) - T? V (L)] ~ L _1 , i.e. there are two correlation length exponents v = 2 and v = 1 that 
govern respectively the averaged/typical loop distribution. 

I. INTRODUCTION 

The stability of pure critical points with respect to weak bond disorder is governed by the Harris criterion : 
near a second order phase transition in dimension d, the bond disorder is irrelevant if the specific heat exponent 
is negative a pure = 2 — dv pure < or equivalently if the correlation length exponent up = v pure > 2/d. On the 
contrary if vp < 2/d, disorder is relevant and drives the system towards a random fixed point characterized by a new 
correlation length exponent satisfying the general bound vp = i>random > 2/d 2] . More recenHy^mportant progresses 
have been made in the understanding of finite size properties of random critical points 0, 0> S El 13 • The main 
outcome of these studies can be summarized as follows. To each disordered sample (i) of size L, one should first 
associate a pseudo-critical temperature T c (z,L), defined for instance in magnetic systems as the temperature where 



the susceptibility is maximum [3J, |fj, |g, |8| ■ The disorder averaged pseudo-critical critical temperature T® V (L) = T c (i, L) 
satisfies 

T? V {L) - T c (oo) - L~ X I VR (1) 

where vp is the correlation length exponent. Eq. Q generalizes the analogous relation for pure systems 

TP ure (L) - T c (oo) - L~ x l Vp (2) 

The nature of the disordered critical point then depends on the width AT C (L) of the distribution of the pseudo-critical 
temperatures T c (i, L) 



AT C (L) = ^Var[T c {i,L)] = ^T*{i,L)-(l 



T c (i,L)) (3) 



When the disorder is irrelevant, the fluctuations of these pseudo-critical temperatures obey the scaling of a central 
limit theorem as in the Harris argument : 

AT C (L) — L~ d/2 for irrelevant disorder (4) 

This behavior was first believed to hold in general 0, , but was later shown to be wrong in the case of random fixed 
points. In this case, it was argued 0,0 that eq. JJJ should be replaced by 



AT C {L) ~ L X I VFL for random critical points (5) 
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i.e. the scaling is the same as the L-dependent shift of the averaged pseudo-critical temperature (Eq. The fact 
that these two temperature scales remain the same is then an essential property of random fixed points that leads to 
the lack of self-averaging for observables at criticality More precisely, for a given observable A, it is convenient 

to define its normalized width as 



x X?(T,L) - (XAT.L)) 2 
Rx(T,L) ee tV - ' v __ v - " (6) 
(*i(T,L)) a 



If £(T) denotes the correlation length, the following behavior of Rx(T, L) is expected 

(i) off criticality, when L>^ (T) , the system can be divided into nearly independent sub-samples and this leads to 
'Strong Self- Averaging' 

Rx {T, L)~ — off criticality for L > £(T) (7) 
Lr 

(ii) in the critical region, when L <C £(T), the system cannot be divided anymore into nearly independent sub- 
samples. In particular at T c (oo) where £ = oo, one can have either 'Weak Self- Averaging' 

Rx(T c {oc),L) ~ L v p for irrelevant disorder (ap < 0) (8) 

or 'No Self- Averaging' 

R x {T c (oo), L) ~ Cst for random critical points (9) 

In this paper, we study from this point of view disordered Poland-Scheraga models with different loop exponents c, 
corresponding to either a pure second order transition with respectively marginal/relevant disorder according to the 
Harris criterion, or to a pure first-order transition. In each case, we numerically compute the histogram of pseudo- 
critical temperatures and study the self-averaging properties at criticality. The paper is organized as follows. In Sec. 
HTl we recall the definition of Poland-Scheraga (PS) models, the critical properties of the pure transitions and the 
disorder relevance. In Sec. 1 1 1 1 1 we describe the observables that we numerically compute for disordered PS models. 
We then present our results for the different loop exponents considered, namely (i) c = 1.5 (marginal disorder) in Sec. 
II VI (ii) c = 1.75 (relevant disorder) in Sec. IY1 fiiil c = 2.15 (case of a first-order transition in the pure case) in Sec I VII . 



II. POLAND-SCHERAGA MODELS : DISORDER RELEVANCE AT PURE CRITICAL POINTS 

A. Critical properties of pure Poland-Scheraga models 

We first consider the adsorption (or wetting) transition of a polymer chain onto an impenetrable substrate. This 
model corresponds to a Poland-Scheraga model with loop exponent c = 1.5. This wetting model in dimension 1 + 1 
is defined as follows. The substrate is located at z = 0. The polymer chain has L monomers, and the position z a of 
monomer (a) satisfies z a > 0, with Z\ — zl = (bound-bound boundary conditions). The partition function of the 
model reads 

z= E ^ H ( 10 ) 

(RW) 

where H = X)q=i £ a ^z a ,o- In equation i|l(J|) . (3 = jrKp is the inverse temperature and the sum runs over all walks 
(RW) such that \z a+ i — z a \ = ±1. This problem can be formulated a la Poland-Scheraga in terms of the loop 
statistics, where each loop of length I has a weight jfj^- We consider here generalized PS models where each loop of 

length I has a weight 4-, with exponent c > 1. 

In the pure case e a = eoj the loop distribution at T c = T c (oo) is 

piro) - ]- c (ii) 

For c > 2, the averaged length < I >= J dl lP^ re (l) is finite, so that the number n(T c ) of contacts with the substrate 
is extensive (n(T c ) ~ L); the transition is therefore first order. For 1 < c < 2, the averaged length < I > diverges, 
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and the Levy sum of n independent variables U drawn from the distribution scales as l\ + ... + l n ~ n 1 /^ *>. As 
a consequence at criticality, the number of contacts n^ ure (T c ) in a sample of length L scales as 

n p L ure (T c ) ~ L c - X (12) 

and the transition is second order. The contact density thus decays as 

pure frp \ 

9 p L ure (T c ) = Ul { c) ~ L c - 2 (13) 
Li 

Its finite-size scaling form 

9 p L ure (T) = L C - 2 Q \(T - T c )lM (14) 
involves the crossover exponent (f>p = = c — 1. 

B. Harris criterion on 2d order transitions 1 < c < 2 

The Harris criterion Q] concerning the stability of pure second order transitions with respect to the addition of 
disorder relies on the sign of the specific heat exponent 

2 C — 3 

ap = 2 — vp = (15) 

c — 1 

An equivalent way to decide whether disorder is relevant consists in a simple power-counting analysis of the disorder 
perturbation exactly at T c : the pure finite-size contact density < S Zit o > plire ~ L c ~ 2 of Eq. (|13(l yields that the 
perturbation due to the presence of a small disorder in the contact energies (4 = eo + Sei scales as 

L 

S *i < S *i,0 >Pure~ L 1 ' 2 x L c ~ 2 = L c ~i (16) 

1=1 

Disorder is thus irrelevant for 1 < c < | and relevant for | < c < 2. The marginal case c — | has been debated for 
a long time [Tol ITU Il2l Il3l 13 ITEl ll^ | and is of special interest since it corresponds to two-dimensional wetting as 
explained above. As a consequence, we have chosen to study in parallel the cases c = 1.5 and c = 1.75. 



C. Disorder relevance on 1st order transitions c > 2 

The effect of disorder on first-order transitions has been discussed for a long time 

El 13 E E3 • In this respect, 

the most frequently studied system 0,123 is the 2D Potts model with q > 4, for which the Aizenman-Wehr theorem 
[l8| states that disorder changes the critical behavior from first order to second order. The recent numerical study of 
this phenomenon j2^, |2^| however shows that there are many subtleties : in particular, the latent heat that vanishes 
for continuous disorder, remains finite for binary disorder. Disorder effects on Potts models have also been studied in 
3D H. 

Since all these studies consider spin systems displaying coexisting domains in the pure case, their conclusions cannot 
be directly applied to the PS-model for the following polymeric reasons [25| : in the PS model, there is no surface 
tension between the localized and delocalized phases, and there exists an diverging correlation length 1/(T C — T) in 
the pure case. So here, the simplest way to discuss the relevance of disorder consists in the simple power-counting 
analysis of the disorder perturbation exactly at T c as in Eq. (|16f) . Here for c > 2 where < S Zii o > pur e is finite, the 
perturbation due to the presence of a small disorder in the contact energies ti = eo + Sei scales as 

L 

53 Se * < 5 *i,° >pure~ L 1/2 (17) 
i=l 

As a consequence, disorder is relevant. 

In this paper, we consider the case c = 2.15 where the pure transition is first order, since this value is of interest for 
DNA denaturation p6l l2?l l28l |29 | . Furthermore, the effect of disorder on this transition has been recently debated 
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III. OBSERVABLES STUDIED IN DISORDERED PS MODELS 



A. Numerical details 

As explained in details in our previous papers the PS models can be numerically studied via the recursion 

relations satisfied by the partition function, with a Fixman-Freire scheme |32j for the entropy of loops. All results 
presented in this paper have been obtained for independent quenched random contact energies (e Q ), distributed 
with a binary distribution. For 1 < c < 2, we have chosen e = or e — Eq with probabilities (1/2,1/2), with 
Eq(c = 1.5) = — 350A and £q(c = 1.75) = — 440A, to obtain critical temperatures in the same temperature range. 
For c = 2.15 of interest in DNA denaturation, we have taken the same values as in our previous study |25j . namely 
e = -355A' or e = -390A with probabilities (1/2, 1/2). 

The data we present correspond to various sizes L with corresponding numbers n s (L) of independent samples. 
Unless otherwise stated, we have considered the following sizes going from L = 2 • 10 3 to L = 2048 • 10 3 , with 
respectively n s (L = 2 • 10 3 ) = 3.84 • 10 6 to n s (L = 2048 • 10 3 ) = 15 • 10 3 . More precisely, we consider 

2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048 (18) 

3840, 1920, 960, 480, 240, 120, 60, 30, 15, 30, 15 (19) 



L 
1000 
n s (L) 
1000 



B. Definition of a sample-dependent pseudo-critical temperature 

In the magnetic systems studied in 0.0]. the pseudo-critical temperature T c (i, L) of the sample i was identified to 
the maximum of the susceptibility. In the PS model, one can not follow the same path and we have tried two different 
definitions : 

Definition of a pseudo- critical temperature from the free- energy 

In the pure PS model with bound-bound boundary conditions, the behavior of the partition function as a function 
of temperature reads 

z pure ,rp rp s , T , T sv P -l 2 L (T C -T)TL 
L i»l/(T-TJ»P l ' 

ZT re (Tc)^^ (20) 

z pure r T>Tc ) ~ 

L V L»1/(T-T C )"-P (T — T C ) 2 L C 

with vp = l/(c— 1). A finite-size pseudo-critical temperature Tf ure ^\L) can then be defined as the temperature 

where the free-energy reaches the extensive delocalized value Fdeioc = — TXln2, i.e. Tc Ure ^\L) is the solution of the 
equation 

ff e (L,r)+TIln2 = (21) 
This definition from the free-energy is very natural, but has the drawback of introducing a logarithmic factor 

/In T \ 1/vp 

T pure U ) {L) _ Tc(oo) „ ( _ \ (22) 

with respect to the purely algebraic factor usually expected (Eq. [21 . This logarithmic factor comes the finite-size 
free-energy value exactly at criticality Ff^ re (L,T c ) = — T C L \n2 + (2 — c)T c InL (Eq. I20|) . In the disordered case, we 

may similarly define a sample-dependent pseudo-critical temperature T c ■ (i,L) as the solution of the equation 

F < f > {L,T)+TL\u2 = Q (23) 

but logarithmic corrections are to be expected at least in the shift (Eq. QJ. Since these logarithmic factors may alter 
the numerical measures of critical exponents in the disordered cases, we have also considered another definition of 
T c (i,L). 

Definition of a pseudo- critical temperature from a sample-replication procedure 
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In the pure PS model with 1 < c < 2, another finite-size critical temperature T< 
temperature where the appropriately rescaled contact densities cross (Eq. I13f) 



■pure(O) 



L 



'[TP ure[e) (L)] = (2L) 2 



rjpure \ 



[T c pMre(e) (2L)] 



(L) may be defined as the 



(24) 



This temperature Tc Ur ' (£) follows the usual algebraic behavior (Eq. |2J, in contrast with the other definition 
Tc Ure (L) that introduces logarithmic corrections (|2~2"|) . 

This definition can be extended to the disordered case via the following strategy first introduced in for the case 
c = 2.15 (with bound-unbound boundary conditions ). In short, one considers a sample (i) of length L, the sample 
(2i) of length 2L, obtained by gluing together two copies of sample i, as well as the sample (Ai) obtained by gluing 
to geth er four copies of sample i : the three contact densities of (i), (2i) and (4i) cross at a temperature (see Fig. 9b 
of |25|) that can be interpreted as the pseudo-critical temperature T c (i, L) of the disordered sample (i) of length L. 

Here we generalize this procedure to the case 1 < c < 2. The basic idea is that in disordered Poland-Scheraga 
models, the loop distribution has exponent c at criticality in the thermodynamic limit as in the pure case l|ll[l . We 
have found numerically this property in our previous studies on the wetting transition with c = 1.5 |l6j | and on the 
DNA denaturation with c = 2.15 j25j- This suggests that the density of contacts scales as L 2 ~ c at criticality as in 

For 1 < c < 2, we have thus tried to define a sample-dependent pseudo-critical temperature T c (i, L) as follows. For 
each sample (i) of length L, we construct the sample (2i) of length 2L by gluing together two copies of sample i, as 
well as the sample (4?) by gluing together four copies of sample i. We then plot the rescaled contact densities 



y T (i,L) = e T (i,L) L 2 



(25) 



with its two analogs y T (2i,2L) = T (2i,2L) (2L) 2 ~ C and y T (Ai,AL) = 9 T (U,AL) (4L) 2 ~ C . Typical results obtained 
respectively for the cases c = 1.5 and c = 1.75 with two disordered samples i = i\, 12 are shown on Fig. ^ : for each 
sample (i), the three rescaled contact densities yr(i,L), yr(2i,2L) and 2/r(4i,4L) cross at a temperature that we 

defined as the pseudo-critical temperature T c (i,L) of the disordered sample i of length L. This crossing of three 
curves validates the present replication procedure to define a proper pseudo-critical temperature. 
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FIG. 1: Illustration of the crossing procedure to obtain the pseudo- critical temperatures T c (i,L) for two samples i = 11,12 (a) 
for the case c=1.5 with L=50000 (b) for the case c = 1.75 with L = 100000. 



Discussion 

Each definition of T c (i, L) has its own advantage and drawback. The definition Tc\i, L), that uses the free-energy, 
is probably the most natural one, but introduces logarithmic corrections, already in the pure PS model. The definition 
T?'(i,L), that uses the sample-replication procedure, may appear more artificial in the disordered case, but it does 
not introduce logarithmic corrections. We have checked numerically that both definitions actually yield very similar 
results for critical exponents and scaling distributions, even though they give different values for a given sample. 
This shows that the conclusions that can be obtained do not depend on the precise definition of the pseudo-critical 
temperature. 
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C. Distribution of pseudo-critical temperatures 

In all cases (c = 1.5, 1.75 and 2.15), and for both definitions of the pseudo-critical temperature (either from the 
free-energy or from the sample-replication procedure) , we numerically obtain that the distributions of pseudo-critical 
temperatures follows the scaling form 

Pl(T c (z, L)) ~ g [x = j (26) 

where the scaling distribution g(x) is simply Gaussian 

g(x) = ^=e~ x ^ 2 (27) 



This means that the only important properties of the pseudo-critical temperatures distribution are the behaviors of 
its average T^ V (L) and width AT C (L) as L varies. 

For the average T® V {L), we have fitted our data with the power-law (Eq^J for the temperature Tc 8 \i, L), and with 
the generalized form involving logarithm (see Eq. 1221 for the pure case) for the temperature T^ [i, L) 

rr W (L) _ Tc(oo) ^-A e {-\ (28) 

/In T \ X I VR 

Tr (/) (i)~r c (oo) ~- A/ (—\ (29) 

The two definitions of T c (i,L) then yield compatible estimates for the exponent vr and the critical temperature 
T c (oo). 

For the variance, we have fitted our data with the power-law (Eq. |SJ for both definitions of pseudo-critical temper- 
atures, and the two definitions of T c (i,L) then yield compatible estimates for the corresponding exponent. 

D. Study of self-averageness in the critical region 

As already mentioned in the Introduction, the issue of self-averageness at random critical points is directly related 
to the scale of the width of the pseudo-critical temperatures : if the width AT C (L) is of the same order of the 

shift T® V (L) — T c (oo), there is no self-averaging at T c (oo). 

We have thus measured for PS-models the L dependence of the ratio Rg(T,L) defined in Eq. © for the contact 
density for temperature T in the critical region 

Re(T,L) = YdMm (30) 
(0i(T, L)) 2 

We have also studied the issue of self-averageness at T c (i, L) via the measure of the ratio 

^ VarMT^ LIL)] 
(0i(T e (i,L),L))» 



We now make contact with the finite size scaling Ansatz of |3j, |6( for an observable X 

X { l\T) = L p Qi ({T-T c {i,L))L 1 ^ (32) 

where sample dependence arises through both T c (i,L) and the scaling function Qi. The ratio Re(T,L) mostly tests 
the relevance of the variance AT C (L), whereas the ratio 1Z C (L) of Eq. I|31|) directly tests the sample dependence 
of Qi(0) for the observable X — 6. We have measured the ratio TZ C (L) for both definitions of T c (i, L). Although 
numerically different, their qualitative behaviour is the same: whenever disorder is relevant, the ratio 1Z C (L) grows 
with L. 

We now present our results for the different loop exponents c. 
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IV. THE CASE OF MARGINAL DISORDER (c = 1.5) 



The mar gina l case c = 1.5 corresponds to two-dimensional wetting and has been the subject of a long-stan ding 
debate [l(J, EH 13 EM 13 • O n the analytical side, efforts have focused on the small disorder limit : Ref [l(| 

finds a marginally irrelevant disorder where the quenched critical properties are the same as in the pure case, up 
to subdominant logarithmic corrections. Other studies have concluded that that the disorder is marginally relevant 
[Til Il2l fl3| . On the numerical side, the same debate on the disorder relevance took place. The numerical studies of 
Ref. |lOj and Ref. have concluded that the critical behavior was indistinguishable from the pure transition. On 
the other hand, the numerical study of ^1 pointed towards a negative specific heat exponent (a < 0), and finally Ref. 
|15| has been interpreted as an essential singularity in the specific heat, that formally corresponds to an exponent 
a = —00. 

It is thus very interesting to study the marginal case c = 1.5 from the point of view of the histogram of pseudo-critical 
temperatures and of self-averaging properties to clarify the situation. 



A. Distribution of pseudo-critical temperatures 



The histograms of both definitions of the pseudo critical temperatures T c (i,L) are shown on Figure E{a) and (b). 
As mentioned previously, both rescaled distributions shown in the Insets of Figure are Gaussian (Eq. 126(1 . 

1.5 i . , . , . 1 1.5 i . , . , . , . 1 



0.5 




FIG. 2: Case c=1.5 : Distribution of pseudo-critical temperatures T c (i, L), together with the rescaling of Eq. 1261 in Inset, (a) 
for definition Tc 6 \i,L) with samples size from L/1000 = 4(0) to 2048 (thick line), (b) for definition T c (/) (i,L) with samples 
size from L/1000 = 2 (dashed line) to 2048 (thick line) 



The averages T^ V (L) behave as 

T2< e \L) ~ 45.15 - 132 f^) 



0.47 



ln^ - 46 



T* vU \L) - 45.25 - 335 f— J (33) 

These values of T c (oo) ~ 45.2 agree with a different determination via the loop statistics in [l^. The mesaured 
exponent vr is very close to the pure exponent vp = 2. However, this does not mean that the disorder is irrelevant, 
since power-law fits of the variances yield 

ATf>(£) -206 ^ 

AT^(L) ^ 188 (^ Z J (34) 

Note that disorder being marginal, extra logarithmic factors could be present in the shift and in the variance. Theo- 
retical predictions on this point are however unavailable, and we have not tried to include these extra factors in the 
above fits. 
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Our conclusion is that the shift and the width have very close exponents, pointing towards a random critical point 
(Eq. |3J) with non self-averaging properties that we now consider. 



B. Non-self-averaging properties 

We show on Figure 01 (a) the behavior of the ratio Rg(T,L) defined in (Eq. 130(1 . Near the critical temperature 
T c (oo) ~ 45.2 obtained in Eq|33| this ratio clearly displays a lack of self averaging, as expected at random critical 
points (Eq. |§J. 

We show on FigureOH (b) the ratios 7Z C (L) of Eq. ((31(1 corresponding to the two definitions of T c (i, L). They both 
grow with L. 



0.12 | . 1 . 1 1 1 . 1 . 1 0.01 




3e+06 



T L 

FIG. 3: Case c=1.5 : (a) Ratio Rg(T,L) for sizes L/1000 = l(O) to 2048 (thick line) (b) The ratios TZ C (L) of (Eq. I3TI . 
corresponding to definitions T^ e \i,L) (thick line) and T^\i,L) (thin line), both grow with L. 



C. Conclusion on the nature of the transition for c = 1.5 



Our numerical results 1(33134(1 indicate that in the presence of disorder, the shift- and width- exponents are close to 
each other and to the pure case exponent — = 0.5. We thus conclude 

^random 2 Impure (*^J 

The disorder is nevertheless relevant in the sense that there is a lack of self-averaging properties at criticality as in 
Eq. @, because both the width AT C (L) and the shift [T c (oo) — T" V (L)] decay with the same exponent 

AT C (L) ~ [T c (oc) - T™(L)] ~ L^' 2 (36) 



V. THE CASE OF RELEVANT DISORDER (c = 1.75) 

A. Distribution of pseudo-critical temperatures 

The histograms of both definitions of the pseudo critical temperatures T c (i, L) are shown on Figure ^a) and (b). 
Both rescaled distributions shown in the Insets of Figure 0] are Gaussian (Eq. 126(1 . 
The respective averaged T^ V {L) behave as 

/ X 0.35 

T^ e \L) ~ 45.2- 155 ( -J 

n t \ 37 

T? V W (L) ~ 45.4 - 140 I J (37) 
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FIG. 4: Case c=1.75 : Distribution of pseudo-critical temperatures T c (i,L), together with the rescaling of Eq. 
(a) for definition T c (e) (i, L) with samples size from L/1000 = l(O) to 128 (thick line), (b) for definition T c (/) (i, L) 
size from L/1000 = 2(Q) to 2048 (thick line) 



126H in Inset, 
with samples 



This shift with the size L shows (Eq. ]QJ that the exponent I/z/r = 0.35 is very different from the pure exponent 
\jvp = c — 1 = 0.75. As expected, the random exponent satisfies the general bound vr>2 Q- For both definitions 
the widths decay with the same exponent, namely 



AT^(L) -134 
AT C (/) (L) -117 



0.38 



0.38 



in agreement with the prediction of Eq. JSJ 



for random critical points. 



AT c (L)~ir(i)-r c (oo) 



Non-self-averaging properties 



(38) 



(39) 



We show on Figure[S](a) the behavior of the ratio i?e(T, L) defined in 13Ufl . It clearly displays a lack of self averaging, 
as expected at random critical points (Eq. . 

We show on Figure[5] (b) the ratios 1Z C (L) of Eq. (|31|l corresponding to the two definitions of T c (i, L). They both 
grow with L. 



Conclusion on the transition for c = 1.75 
For relevant disorder c = 1.75, our conclusion is that the width AT C (L) and the shift \T c (po) — T°"(L)] decay with 



2.7 > 2 > Vpure) and there is no self-averaging at criticality, in 



the same new exponent L ^/"random ( w here ^random 

agreement with the general predictions J5j and © for random critical points. 

VI. DISORDER EFFECTS ON A FIRST ORDER TRANSITION (c = 2.15) 

The case c = 2.15 is of special interest in the context of DNA denaturation, and the nature of the transition in the 
presence of disorder has been under debate recently [2j| |3(1 . The numerical studies |2^, l3fl| with bound-unbound 
boundary conditions have found crossings of the energy density (Fig. 4 of [3(| and Fig. 7b of [25|) and of the contact 
density (Fig. 6 (b) of [25|) for various sizes L. These results point towards a finite energy density and a finite contact 
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density at criticality. However, a recent probabilistic analysis |3jJ of the disordered PS model with c > 2 has concluded 
that the second derivative of the free-energy remains bounded when approaching the critical point from the localized 
phase, and that the order parameter vanishes continuously. It is thus interesting to reconsider the problem from the 
point of view of the distribution of pseudo-critical temperatures to clarify the nature of the transition. 



A. Distribution of pseudo-critical temperatures 

The histograms of both definitions of the pseudo critical temperatures T c (i,L) are shown on Figure^a) and (b). 
Both rescalcd distributions shown in the Insets of Figure H3 are Gaussian (Eq. I26JI . 





FIG. 6: Case c=2.15 : Distribution of pseudo-critical temperatures T c (i,L), together with the rescaling of Eq. I|26|l in Inset, 
(a) for definition Ti [i, L) with samples size from L/1000 = l(O) to 128 (thick line) (b) for definition T c (/) (i,L) with samples 
size from L/1000 = 2 (dashed line) to 512 (thick line) 



The averages T^ V (L) behave as 



81.36- 



1254 



L 



81.38- 357 



InL 



(40) 



The two values of T c (oo) are close and agree with the value obtained in Ref. [25| with a different method. On one 
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hand, according to Eq. ^ this behavior of the shift seems to indicate that the exponent is unchanged with respect to 
the pure case Vr = 1 = Vp. On the other hand, the width is found to be much bigger than the shift since 



AT c (e) (L) ~ 23.8 (\ 



L 



0.4!) 



0.49 



ATP(L) ^23.2(j-j (41) 
This behavior corresponds to the Central Limit estimation L" 1 / 2 , as in the Harris argument. 

The fact that the shift and the width exhibit different exponents actually means that two diverging correlation 
lengths coexist in the presence of disorder 



£,shift{T) 



1 



(T c - T) 



S,var{T) ~ .y, ^ ^ 2 (42) 

The presence of two different correlation length exponents was already found in other models, in particular in the 
random transverse field Ising chain |33j , where the exponent v = 2 governs the decay of averaged correlation, whereas 
v = 1 governs the decay of typical correlations. In Sec. VII A of Ref. [3i|, a scaling analysis of disorder effects on 
first-order phase transitions in dimension d suggests that the presence of two different correlation length exponents 
v = 2/d and v = 1/d should be generic in these systems : the exponent v = l/dis expected to describe the rounding 
of the transition in a typical sample, whereas v = 2/d describes the rounding of the transition of the distribution of 
samples. We will discuss this issue in more details in Section IVl Dl after the description of our numerical results on 
non-self-averaging properties. 

B. Non-self-averaging properties 

We show on Figurc[7|(a) the behavior of the ratio Re(T, L) defined in 13U|) . It clearly displays a lack of self averaging 
at criticality. 

We show on Figure [7| (b) the ratios TZ C (L) of Eq. (|31fl corresponding to the two definitions of T c (i, L). They both 
grow with L. 




3e+06 



FIG. 7: Case c=2.15 : (a) Ratio R e (T,L) for sizes i/1000 = l(O) to 2048 (thick line) (b) The ratios TZ C {L) of (Eq. I3TI . 
corresponding to definitions Tc 6 \i, L) (thick line) and T^\i,L) (thin line), both grow with L. 
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C. Finiteness of the contact density at criticality 

The dominance of the variance (|41|l over the shift l|40|) indicate that asymptotically for large L, half of the samples 
(i, L) are still localized at T c (oo), whereas the other half is already delocalized. This suggests that the contact density 
is finite at criticality, as we have numerically found in our previous study [25J. However, since Ref. [3l| states that 
the order parameter vanishes continuously at the transition, we have performed more detailed calculations. Figure [S] 
(which is more detailed than Figure 6 (b) of our previous study 25] ) shows the averaged contact density 6l(T) as a 
function of the temperature T for several sizes : the results cross regularly without rescaling. This points towards a 
finite contact density at criticality. 




81.1 81.2 81.3 81.4 81.5 81.6 

T 



FIG. 8: Case c=2.15 : crossing of the averaged contact density 6l(T) for sizes from L/1000 = l(O) to 2048 (thick line) . 



D. Discussion on the nature of the transition for c = 2.15 in the presence of disorder 

Our present study shows that, in the presence of disorder, the transition for c = 2.15 is an unconventional random 
critical point with two different correlation length exponents v = 2 and v = 1 l|42l) . This is in contrast with usual 
random critical points, arising from second order transitions with relevant disorder, where the same exponent is 
expected to gov ern the width and the shift (Eq. 0), but this is reminiscent of what happens at strong disorder 
fixed points |33l l34| . The question is now which correlation exponent appears in a given observable. In the random 
transverse field Ising chain where many exacts results are known for exponents and scaling distribution functions |33| , 
it is well understood how the two exponents v = 2 and v = 1 govern respectively the averaged/typical correlations. 
Here in the disordered PS model, the analog of the correlation function is the loop distribution. To simplify the 
discussion, let us more specifically consider the probability of an end-to-end loop of length L in sample (i) of length 

L, which is directly related to the partition function Z^'(T) of sample (i) 

pU(L,T) = ^-J— (43) 

Introducing for each sample (i) the difference between the free-energy density F ( /\T)/L = -T\nZ^{L,T)/L and 
the delocalized value fdeioc = —Tin 2 

/(-)(L,T)^- rlnZW ^ r ) + rin2 (44) 
Li 

one gets 

lnP£° (L,T) = — c hxL + L(3f^ (L, T) (45) 

The self-averaging property of the free energy means that f^(T) converge for large L to a non-random value f(T) 
for any sample (i) with probability one 

ifm — f(T) (46) 
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where f(T) is the free-energy difference between the localized phase and the delocalized phase : f(T < T c ) < and 
f(T > T c ) = 0. This translates immediately into the corresponding statement l|45[) for the logarithm of end-to-end 
loop probability 

^^l^fm (47) 

for any sample (i) with probability one. Since the typical correlation length £(T) is usually defined as the decay rate 
of the logarithm of the correlation, we obtain here that it is simply given by the inverse of the free-energy f(T) 

1 s . )im (^n) = . mT) m 



£(t) l^oo y L 

Here for c = 2.15, we expect from the discussion of the previous Section fVI CI that the contact density is finite at 
criticality : this implies that the energy is also finite at criticality, and thus we are led to the conclusion that the 
free-energy vanishes linearly 

f(T) — » (T c -T) (49) 

T->T C ~ 

The typical correlation length involves the exponent v = 1 

i(T) — ► (T c - T)-" with v = 1 (50) 

T—*T~ 

Let us now consider the decay of the averaged end-to-end loop distribution that defines an a priori different 
correlation length £(T) 



ln(>f (L,T) 



L L^oo £(T) 



(51) 



This correlation length £(T) determines the divergence of high moments of the averaged loop distribution. At a given 
temperature T < T c , these moments will actually be dominated by the rare samples of length L which are already 
delocalized at T, i.e. the samples having T c (i,L) < T. Since our numerical results indicate that the distribution of 
the pseudo-critical temperature T c (i, L) is a Gaussian with mean and width given respectively by Eqs (|40|) and H41JI. 
we obtain that the fraction of delocalized samples presents the following exponential decay in L 

Vmh[T c (i, L) < T] ~ e -( T c°°- T ) 2L (52) 

This measure of the rare delocalized samples will govern the decay of the averaged loop distribution, and the correlation 
length defined in l|51|) thus involves the exponent v = 2 

£(T) — ► (T c - T)- v with v = 2 (53) 

T— T~ 

in contrast with the typical correlation length (|50|l . 

To better understand the emergence of two different correlation lengths, we have numerically measured the distri- 
bution over the samples (£) of the free-energy /W(L, T) defined in Eq. H44|) . We obtain that for T < T c 

/»(D=/(D + ^ + ^ (54) 

where ax is temperature dependent and Ui is a Gaussian random variable of zero mean and of variance 1 

G{u) = -±=e-^ (55) 

V Z7T 

The averaged end-to-end loop distribution then reads l|45|) 

I 1 /"-t-OO -1 rt 2 rr 2 

pf(L,T) = J- e W ( «(£,r) = J_ e W(D / du G(u)e ^<r r « = i e W(T)wV ( 56 ) 
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The difference between the correlation length £(T) l|51|) and the typical correlation length (|48|l is due to the variance 

2 

£(T) £(D 2 1 j 

In particular, to obtain the scaling ~ (T c — T) 2 different from — 1_ ~ (T c — T), the variance term in a\ has to 

cancel exactly the leading order in (T c — T) on the left hand-side. 

So the picture that emerges of the present analysis is very reminiscent of what happens at strong disorder fixed 
points H3, 0| : the exponents v — 1 and v = 2 govern respectively the decay of typical/ averaged loop distribution. 
Our conclusion is thus that the exponent v = 1 governs the free-energy l|49|l that corresponds to a Lyapunov exponent, 
i.e. it describes the critical behavior of any typical sample, whereas the exponent v = 2 = 2/d is the finite-size scaling 
exponent of Chayes et aZ|3 and is related to the variance of the distribution of pseudo-critical temperatures. Figure 
8 of our previous paper |25j may be now interpreted as follows : for each sample, the critical region has a width of 
order whereas the contact density averaged over the samples decay on a much wider scale that represents 

the sample-to-sample fluctuations of the pseudo-critical temperatures T c (i, L). 



VII. SUMMARY AND DISCUSSION 

In this paper, we have used the recent progresses in the theory of finite-size scaling in disordered systems 0, 0, IE 
0- 0- 0] t o study the role of disorder in Poland-Scheraga models. We have obtained that the numerical measure of 
the distribution of pseudo-critical temperature T c (i,L) over the samples is a very powerful tool to elucidate the true 
nature of random critical points. The comparison between the averaged shift T^ V (L) — T c (oo) and the width AT C (L) 
clarifies the role of disorder, and allows one to understand the non self averaging properties of various observables at 
T c (oo) whenever disorder is relevant. 

For c = 1.75 corresponding to a second order transition with relevant disorder, we have obtained that both the width 
AT C (L) and the shift [T c (oo) — T£ V (L)] decay as L~ 1 ^ random with the same new exponent v r andom ~ 2.7 > 2 > v vure . 

For c = 1.5 corresponding to a second order transition with marginal disorder, we have obtained that the exponent is 
unchanged v — 2 with respect to the pure case, but that disorder is nevertheless relevant : the width AT C (L) ~ L^ 1 / 2 
is of the same order of the shift T^ V (L) — T c (oo) ~ L -1 / 2 and this means that thermodynamic observables (apart 
from the free-energy) remain distributed at criticality. We have checked it numerically for the contact density. 

For c = 2.15 corresponding to a first order transition in the pure case, we have obtained that the width AT C (L) ~ 
L~ x / 2 dominates over the shift T® V (L) — T c (oo) ~ L _1 . The presence of two correlation length exponents v = 2 and 
v = 1 is reminiscent of what happens at strong disorder fixed points |33l l34j , and we have explained how these two 
exponents appear respectively in the typical and averaged loop distributions. 
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